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Abstract 

We review probabilistic approaches to the Gross-Pitaevskii theory de- 
scribing interacting dilute systems of particles. The main achievement 
are large deviations principles for the mean occupation measure of a large 
system of interacting Brownian motions in a trapping potential. The cor- 
responding rate functions are given as variational problems whose solution 
provide effective descriptions of the infinite system. 



1 Introduction 



The phenomenon kno wn as Bose-Eins tein condensation (hereafter abbreviated 
BE C) was predi cted bv lEinsteinI (|l925l ) on the basis of ideas of the Indian physi- 
cist Bosd ( 1924 ) concerning statistical description of the quanta of light: In a 
system of particles obeying Bose statistics and whose total number is conserved, 
there should be a temperature below which a finite fraction of all the particles 
"condense" into the same one-particle state. Einstein's original prediction was 
for a non-interacting gas of particles. The predicted phase transition is asso- 
ciated with the condensation of atoms in the state of lowest energy and is the 
consequence of quantum statistical effects. 

For a long time these predictions were considered as a curiosity of non- 
interacting gases and had no practical impact. After the obser vation of super - 
fluidity in liquid '*IIe below the A temperature (2.17 K) was made, London ( 19381 ) 
suggested that, despite the strong interatomic interactions, BEC indeed occurs 
in this system and is responsible for the superfluidity properties. This sugges- 
tion has stood the test of time and is the basis of our modern understanding of 
the properties of the superfluid phase. 

The first self-consistent theory of super-fluids was developed by iLandau 



(|l941r ) in terms of the spectrum of elementary excitations of the fluid. In 1947 
Bogoliubov developed the first microscopic theory of interacting Bose gases, 
based on the concept of Bose-Einstein condensation. This initiated several 
theoretical studies; a recent acco unt on the state of the art can be found in 
Adams. S. and Bru. J.-B. I (l2004a|lbl'l and on its contribution to superfluidity the- 
orvin lAdams. S. and Bru. J.-B.I (|2004cl ). After lLandau. L.D. and Lifshitz. E.M 



(1951) had appeared, Penrose! (1951) and Onsager. L. and Penrose. O. ( 1956 ) 



introduced the concept of the non-diagonal long-range order and discussed its re- 
lationship with BEC. An important development in the f ield took p l ace w ith the 
prediction of quantised vortices bv lOnsager ( 1949| ) and Fevnman ( 1955 ). The 
experimental studies on dilute atomic gases were developed much later, starting 
from the 1970s, benefiting from the new techniques developed in atomic physics 
based on magnetic and optical trapping, and advanced cooling mechanisms. 

In 1995, the first experimental realisations of BEC were achieved in a system 
that is as different as possible from ^He, namely, in dil ute atomic alkali gases 



trappe d by magnetic fields. These re alisations are due tolAnderson. M.H. et al 



(|l995l ). iBradlev. C.C. et al.l (|l995[) . iDavis. K.B. et all (ll995|), after appropri- 
ate cooling methods had been developed. For this remarkable achievement, 
the Nobel prize in physics 2001 was awarded to E.A. Cornell, W. Ketterle 
and C.E. Wieman. Over the last few years these systems have been the sub- 
ject of an explosion of research, both experimental and theoretical. A com- 
prehensive account on Bose-E i nstein condensation is the recent monograph 



Pitaevskii. L. and Stringari. S.I (|2003[) 



Perhaps the most fascinating aspect of BEC is best illustrated by the cover 
of Science magazine of December 22, 1995, in which the Bose condensate is 
declared as the "molecule of the year" . The Bose condensate is pictured as a 
platoon of soldiers marching in lookstep: every atom in the condensate must 
behave in exactly the same way. One of the most striking consequences is that 
effects, which are so small that they are practically invisible at the level of a 
single atom, are spectacularly amplified. 

Mot i vated by the experiment a l success, in a series of p a pers Lieb. E.H. et al.l 

( 2000al ). lLieb. E.H. et al.l (|2000bh . lLieb. E.H. et a"n(|200l[ ). lLieb. E.H. and Seiringer. R 
( 2002f ) obtained a mathematical foundation of Bose-Einstein condensation at 
zero temperature. The mathematical formulation of the A'^-particle Boson sys- 
tem is in terms of an TV-particle Hamilton operator, Hn, whose ground states 
describe the Bosons under the influence of a trap potential and a pair potential, 
see Section [51 Lieb et al. rigorously proved that the ground state energy per 
particle of Hn (after proper rescaling of the pair potential) converges towards 
the energy of the well-known Gross-Pitaevskii functional. The ground state is 
approximated by the iV-fold product of the Gross-Pitaevskii minimiser mulit- 
plied by a correlated term involving the solution of the associated scattering 
equation. Moreover, they also showed the convergence of the reduced density 
matrix, which implies the Bose-Einstein condensation. As had been generally 
predicted, the scattering length of the pair interaction potential plays a key role 
in this description. 

These rigorous results are only for zero temperature, whereas the experi- 
ments show BEC at very low, but positive temperature. The mathematical 
understanding of BEC at positive temperature is rather incomplete yet. Its 
analysis represents an important challenging and ambitious research area in 
the field of many-particle systems. Thermodynamic equilibrium states are de- 
scribed by traces of e"'^^" , where /3 G (0, oo) is the inverse temperature and Hn 
is the A^-particle Hamilton operator. Via the Feynman-Kac formula (see e.g. 
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FevnmaiJ (Il953l) and lGinibrd dHz^)), these traces are expressed as exponential 



expectations of N interacting Brownian motions with time horizon [0,/3]. This 
opens up the possibihty to use probabilistic approaches for the study of these 
traces, in particular stochastic analysis and the theory of large deviations. 

In this review we present our probabilistic approaches to dilute systems 
of interacting many-particle systems at positive temperature using the Gross- 
Pitaevskii approximation. Using the the theory of large deviations, we char- 
acterise the large- iV and the large-/? behaviour of various exponential expec- 
tations of N interacting Brownian motions with time horizon [0, 0\ in terms 
of variants of the Gross-Pitaevskii variational formula. In particular we intro- 
duce and analyse a new model, which we call the Hartree model, whose ground 
states are the ground product states of the Hamilton operator Hn- Their large- 
behaviour is characterised in terms of the Gross-Pitaevskii formula, with 
the scattering length replaced by the integral of the pair interaction potential. 
This nice assertion is complemented by an analogous result for positive tem- 
perature. Our programme started with lAdams. S. et al. (l2006allbl'l. which we 
su mma rise here. 
(20073). 



Further aspects are con si dered in lAdams. S. and Dorlas. T 



Adams. S. and Konig. W.I (|2007f) . lAdams. S. and Dorlas. T l ( 2007bf) 



and I Adams' (200 7a|). Under cu rren t development ar e lAdamd (|20073) . lAdamsl 
( 2007d) . , Adams. S et al.l (|2007l ) and l Adams. S. et all (l2007|) in which non-dilute 
systems are studied. 

We give a brief introduction to the physics of dilute quantum gases and 
their mathematical treatment at zero temperature in Section [2l In particular we 
introduce the Gross-Pitaevskii formula and the scattering length and describe 
the results by Lieb et al. and our results of the ground product state. Our 
probabilistic models are introduced in Section [3l Section [4] is devoted to our 
large deviations results and the variational analysis. 



2 Dilute Quantum Gases 

We introduce the modelling of the Gross-Pitaevskii theory which will be the 
starting point for our probabilistic models in Section [3) Let us comment briefly 
on some issues of the 1995 experiments as these are the motivation for the 
renewed interest in the Gross-Pitaevskii theory and its analytical proof by Lieb 
et al. 

The experimental systems are collections of individual neutral alkali-gas 
atoms (e.g., ^Li, ^qk, ^TRb, 23Na, ^Li and **5Rb, *^7Rb,i33Cs,i''4Yb, ^^Rh2, 
and ^Li2)), with total number N ranging from a few hundreds up to ^ 10^°, 
confined by magnetic and/or optical means to a relatively small region of space. 
Their densities range from ~ lO^^cm^^ to 5 x lO^^cm^^, and their tempera- 
tures are typically in the range of a few tenths of nK up to ~ 5/j,K. 

In a typical system, we are faced with several length scales. One of them 
is the two-body interaction energy Y? jmo?, where m is the reduced mass of 
the two particles, h, is Heisenberg's constant and a is the scattering length 
(see Section [2^1 below), expressing the strength of the interatomic interaction. 



3 



A second one is the mean interparticle spacing r-mt, and a third one is the 
oscillator frequency Oosc of the confining trap potential. Note that the first scale 
does not depend on the trap geometry, whereas the oscillator frequency aosc, 
the mean interparticle spacing, the transition temperature Tc and the mean- 
field energy Uq (to be specified later) do depend on the shape of the confining 
potential. Introduce the "healing length" ^ = {2mnUoh)~^/'^ and the de Broglie 
wavelength Adb- Note that Oosc is the zero-point spread of the ground-state 
wave function of a free particle in the trap. The relations between these scales 
are as follows. 

a < Hnt ^ Adb < < Oosc- 

Typical values are a 50 A, Tint ^ 2000 A, ^ ^ 4000 A, flosc ^ 1m- If one 
compares these numbers with those of liquid helium, one sees that the dilute 
gas condition a <^ Tint, which is characteristic for the BEC of alkali gases, is 
very far from satisfied for liquid helium. As a consequence, liquid helium is 
a much more strongly interacting system than BEC gases, by many orders of 
magnitude. 

We now turn to a mathematical modelling and introduce the potentials 
and the scattering length in Section 12.11 and the Gross-Pitaevskii theory in 
Section!^ 



2.1 Potentials and Scattering Length 

Our two fundamental ingredients are a trap potential, W, and a pair-interaction 
potential, v. Wc restrict ourselves to dimensions d € {2,3}. Our assumptions 
on W are the following. 

W: M."^ ^ [0, oo] is measurable and locally integrable on {W < oo} with 
lim inf W{x) — oo. 

-R— >CXD \x\>R 

In order to avoid trivialities, we assume that {W < oo} is either equal to K'' or 
is a bounded connected open set containing the origin. 

Our assumptions on v arc the following. By Br{x) we denote the open ball 
with radius r around x G K'^. 

V : [0, oo) — + M U {-t-oo} is measurable and bounded from below, 
a :— sup{r > 0: v{r) — oo} G [0, oo), 'f^l[?),oo) is bounded Vt^ > a. 

Note that we also admit v{a) — +oo. We are mainly interested in the case 
where v has a singularity, i.e., either a > 0, or a = and limr^o ''^('') = oo- 
Exa mples include a lso super-stable potentials and potentials of Lennard- Jones 
type ( Ruelld ( 1969^ )). According to integrability properties near the origin, we 



distinguish two different classes as follows. We call the interaction potential 
V a soft-core potential if a — and Jg^^f^^v{\x\) dx < +oo. Otherwise (i.e., if 

a > 0, or if a = and /g^j-o) ''^(Is;!) dx = +oo), we call the interaction potential 
a hard-core potential. 
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We shall need the following dAf-dimcnsional versions of the trap and the 
interaction potential: 



N 

W{x) ^'^W{xi) and = ^ v{\xi - Xj\), 

8=1 l<i<j<N 

where x = {xi, . . . , xjy) G . 

Let us introduce the scattering length of the pa ir potential, v, and its most 



impor tant properties. For a detailed overview, see lLieb. E.H. and Yngvason. J 



(|2001l ). First we turn to d > 3. Let u: [0, oo) — > [0,oo) be a solution of the 
scattering equation, 

u" = ^uv on (0, oo), u{0) = 0. (3) 
Then the scattering length a(v) G [0, oo], of u is defined as 

a{v) = hm r j-A- . (4) 

1 — ►oo L u [r) i 

If v{0) > 0, then a{v) > 0, and if v{r)r'^~^ dr < oo, then a{v) < oo. In the 
pure hard-core case, i.e., v = ooII^q we have a{v) — a. It is easily seen from 
the definition that the scattering length of the rescaled potential ^~^v{-£,^^) is 
equal to £,a{v), for any ^ > 0. 

There is some ambiguity of the choice of u in ([3]); positive multiples of u 
are also solutions, but the factor drops out in We like to normalise u by 
requiring that limf(^aou'{R) = 1. It is easily seen that (where Ud denotes the 
area of the unit sphere in R''), 

I- — dx — LUd / v{r)u{r)r dr — 2ujd / u"{r)rdr 

m Jq Jo 

= 2ujd lim (u{r)r — / u (r) dr\ 

-R^ooV Jq J 

= 2ijJd lim (u'{R)R ~ u{R)) = 2ujda{v). 

As a consequence, in dimension d = 3, we have a(v) < a{v). Indeed, u 
is a nonnegative convex function whose slope is always below one because of 
linii^^oo u'{R) = 1. By m(0) = 0, we have that u{r) < r = r'^~'^ for any r > 0. 
With the help of ^ we therefore get 87rQ;(t;) = 2uoda{v) < ^^iv{\x\)dx = 
87rS(ti). 

In d = 2, the definition of the scattering length is slightly different. We treat 
first the case that supp(w) C [0, i?*] for some i?* > and consider, for some 
i? > i?*, the solution u: [0, i?] [0, oo) of the scattering equation 

u" = ^uv on[0,i?], = l,'u(0) = 0. 
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Then u{r) = ^^S sfi^ ^ ^'^^ i?, < r < i? for some a{v) > 0, which is by 

definition the scattering length of v in the case that supp(w) C [0, i?,]. Note 
that a(v) does not depend on R. Hence, 

log r-u(r) log i? T, ^ ^ T} 

loga w = — , R^<r<R. 

1 — u[R) 

For general v (i.e., not necessarily having finite support), v is approximated by 
compactly supported potentials, and the scattering length of v is put equal to 
the limit of the scattering lengths of the approximations. 

The dilute gas condition ensures that the scattering length is a satisfactory 
measure of the interaction strength. This approximation neglects any higher 
energy scattering processes. We finally discuss briefly the effects of the atom- 
atom scattering on the properties of the many-body alkali-gas system. The 
fundamental result is that under some conditions the true interaction potential 
V of two atoms of reduced mass m may be replaced by a delta function of 
strength 2iTt^a/m. The effective interaction is 

v^s{x) ^ S{x), X G R"^. 

m 

This motivates to scale the potential in such a way that it approximates the 
delta function in the large A^-limit. This will be done in the so-called Gross- 
Pitaevskii scaling in Subsection 12. 2[ which is a particular approximation of the 
delta function. 



2.2 The Gross-Pitaevskii approximation 

The simplest possible approximation for the wave function of a many-body 
system is a (correctly symmetrised) product of single-particle wave functions, 
i.e., the Hartree-Fock ansatz. In the case of a BEC system at temperature 
T = 0, this approximation usually leads to the Gross-Pitaevskii approximation. 
Basically the Gross-Pitaevskii approximation suggest to replace the evolution 
(time-dependent or stationary) of the many-body wave functions, governed by 
a system of Schrod i nger e quations, by a q ne-particlc non-linear Schrodinger 
equation (see GrossI ( 196ll) . IPitaevskiil ( 196l[ )): 



idt'^{x,t) = - + M/ + 47ra|^'(a;,i)|^j*(a;,t), x eW.te R+. 

In the stationary case the Gross-Pitaevskii theory gives an approximation for the 
quantum mechanical ground state for many particles (i.e., in the limit N —>■ oo) 
as a variational problem for a single particle in an effective potential. Hence we 
first summarise some ground state properties for finitely many particles. 
The ground-state energy per particle of the A^-particle Hamilton operator 
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is given by 

Here H^{R'^) = {/ G L^{R'^) : V/ G L2(R'^)} is the usual Sobolev space, and 
V is the distributional gradient. It is standard to proof that there is a unique, 
continuously diffcrentiablc, minimiser h^, £ H^{M.'^) on the right hand side of 
([6]), and that it satisfies the variational equation 

AK = + o/i* - NxnK- 

Now we turn to the above mentioned product ansatz. Introduce the ground 
product state energy of Ti^r, that is, 

= -T-r inf (/ii «)■•■«) ft. AT, Hjv/ii «)■•■«) /ijv). (7) 

N hi,...,hNeH^(M.'i): \\hi\\'i=iyi^ ' 

The replacement of the ground state energy, xn, by t he ground product state en- 

ergy. y '^,? ' , is known as the Hartree-Fock approach (see lPickhoff. W.H. and Van Neck. D 
(j2005l )). Sometimes, the formula in ([7]) is called the Hartree formula. Obviously, 

Xn ^ Xn- 

We can also write 

N 



x7 



iV hi hjvEff (« ^ ^ J — J 



l\hi |h = l Vi 



l<i<j<N 



where V denotes the integral operator with kernel w o | ■ | , either defined for func- 
tions by V^/(a;) = /jjj v{\x-y\)f{y) dy or for measures by V^^(x) = J^a /i(dy) v{\x' 
y\). The main assertio ns on the formula in ([7]) and its minimisers are summarised 
as follows (sec Adams, S. et al . (2006a)). 

Lemma 2.1 (Ground product states of ?tjv) Fix G N. 

(i) There exists at least one minimiser (fti, . . . , ft^r) of the right hand side 
in the formula for xJv ' ■ '^he set of minimisers is compact and invariant 
under permutation of the functions hi, . . . , h^ ■ 

(a) Any minimiser {hi, . . . , /iat) satisfies the system of differential equations 

Ah^ = -\^h^ + Wh^ + h^Y^ Vk^^, i=\,...,N, 

With X, = W^KWl + {W,hf) +Y.J^^{hhyh'^)■ Furthermore, \\h,\\^ < 
Cd{\^ - (iV - l)infw)''/4 for any i G {1, . . . , N}, where Cd > depends 
on the dimension d only. 
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(in) Let v be soft-core, assume that d E {2,3}, and let (hi, . . . , hj^) be any 
minimiser. Assume that w|(o,?7) ^ for some 77 > 0. In d = Z, furthermore 
assume that 

/ |''^(l2/l)|^~''*d2/ < 00, for some 5 > Q. 

Then every hi is positive everywhere in M'' and continuously differentiahle, 
and all first partial derivatives are a-Holder continuous for any a < 1. 

(iv) Let V be hard-core, assume that d G {2,3}, and let (hi, . . . , h^) he any 
minimiser. Then every hi is continuously differentiahle in the interior of 
its support, and all first partial derivatives are a-Holder continuous for 
any a < 1 . 

Remark 2.2 (i) Unlike for the ground states ofTi.jsf in there is no con- 
vexity argument available for the formula in ([7]). This is due to the fact 
that a convex combination of tensor-products of functions is not tensor- 
product in general, and hence the domain of the infimum in ([7]) is not a 
convex subset of H^{W^^). However, for /12, • ■ ■ , ft-TV fixed, the minimisa- 
tion over hi enjoys the analogous convexity properties on H^{W^) as the 
minimisation in 

(a) If v is hard-core, it is easy to see that the distances between the supports 
of hi, ... , hj^ have to be no smaller than a (see ^) in order to make the 
value of {hi®- ■ -^hiq ,T-Li^hi®- ■ -^hiq) finite. The potential J^j^i^^^ 
equal to 00 in the a-neighbourhood of the union of the supports of hj with 
j ^ i, and hi is equal to zero there (we regard Q ■ 00 as Q). In particular, 
minimisers of ([7]) are not of the form {h, . . . ,h). In the soft-core case, 
this statement is not obvious at all. A partial result on this question in 
d — i will be a by-product of Section \2.2\ below. 

O 

We study now our main variational formulas, xn and x^v'' their minimisers 
in the limit for diverging number TV of particles. In particular, we point out 
some significant differences between xn and its product state version Xn' 
the soft-core and the hard-core case, respectively. 

First we report on recent results by Lieb, Seiringer and Yngvason on the 
large- iV behaviour of xn- Let the pair functional v be as in ([2|) and assume 
additionally that v >Q and v{0) > 0. 

We shall replace v by the rescahng VAr(-) — £,'^'^v{ ■ C^^), for some appropriate 
^AT tending to zero sufficiently fast. This will provide the dilute gas condition 
needed. Hence, the reach of the repulsion is of order ^at, and its strength of 
order Furthermore, the scattering length of u, a{v), is rescaled such that 
a(t^Af) — cx.{v)I3m- If (^N i sufficiently fast, this rescahng makes the system 
dilute, in the sense that a{vN) < N-^/"^. This means that the interparticle 
distance is much bigger than the range of the interaction potential strength. 
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More precisely, the decay of 13n will be chosen in such a way that the pair- 
interaction has the same order as the kinetic term. 

The mathematical description of the large- TV behaviour of xn in this scal- 
ing, and hence the theoretical foundation of the above mentioned physical 
experiments, ha s been successfully accomplished in a recen t series of papers 
Lieb E H. et all f2 000al. Licb, E. H and Yngvason. J.I (|200l[) . iLieb. E.H. et a\ 



(|200lh . iLieb. E.H. and Seiringer. R.I (|2002[ ). It turned out that the well-known 
Gross-Pitaevskii formula adequately describes the limit of the groun d state s and 
its energy. This variational formu la was first intro duced in iGrosd (Il96lh and 
Gross! ( 1963 ) and independently in iPitaevskii ( 1961 ) for the study of superfluid 
Helium. After its importance for the description of Bose-Einstein condensation 
of dilute gases in magnetic tra ps was realised in 1995, the interest in this for- 
mula consi derably increased; seelDalfovo. F. et al. I (|l999l) for a summary and the 
monograph Pitaevskii. L. and Stringari. S.I ( 2003 ) for a comprehensive account 



on Bose-Einstein condensation. 

The Gross-Pitaevskii formula has a parameter a > and is defined as follows: 



y(gp) 



inf {\\Vm + {W,^'')+AnaU\\\). 



It is known Lieb. E.H. et al.l ( 2000a[ ) that xL*^^' possesses a unique minimiser 
4>''a^\ which is positive and continuously differentiable with Holder continuous 
derivatives of order one. 

Since w(0) > 0, its scattering length a{v) is positive. The condition 



v{r)r dr < oo 



a+l 



implies that Q!(u) < oo. Furthermore, note that the rescaled potential ^ '^^{■^ ^) 
has scattering length £,a(v) for any ^ > 0. 



Theore m 2.3 (Large-A/' asymptot i c of y jy in d £ |2, 3} ) iLieb. E.H. et al 
t2000a) , Lieb. E.H. and YngvasonTj] (|200ll) , iLieb. E.H. etail (l200ll) ]. A 

that d e {2, 3}, that v>0 with v{0) > 0, and v{r)r'^~'^ dr < oo. Repl 
byvN{) = S,N^v{ ■ Civ^) = ^/N ind — 3 and — a{v)~ 

in d = 2. Let hjM e H^IW^^) be the unique minimiser on the right hand side of 



ssume 
ace V 



([6|). and define (f>ff G as the normalised first marginal o/hj^, 



4>\,{x) = / h\{x,X2, ■ ■ ■ TXN)dx2- ■ ■ dxN, x^W^. 

Then we have 

^lirn^Xw = Xa(v] <^Af i'Pai'v))'^ L'^ {M.'^)-sense. 

In particular, the proofs show that the ground state, /ijv, approaches, for 
large N, the function 

{xu...,xn)^Y[ ( iu(gp)ii' /(mindxi -XjI: j < i})\ 
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where /(r) — u{r)/r and u is the solution of the scattering equation ([3]). In 
order to obtain the Gross-Pitaevskii formula as the hmit of XAf also in d = 2, the 
rescaling of v in Theoreni l2.3l has to be chosen in such a way that the repulsion 
strength is the inverse square of the repulsion reach and such that this reach 
decays exponentially, which is rather unphysical. 

Ther e is an analogue of Theor em 12.31 for the Hartree model in the soft-core 
case, see I Adams. S. et aL ( 2006a ). It turns out that the ground product state 



energy x^' also converges towards the Gross-Pitaevskii formula. However, in 
c? = 2, it turns out that the potential v has to be rescaled differently. Further- 
more, in d G {2, 3}, the scattering length a{v) is replaced by the number 

Theorem 2.4 (Large- iV asymptotic of X^n' i soft-core case) Letd E {2,3}. 
Assume that v is a soft-core pair potential with v > and v{0) > and 
a{v) < oo. In dimension d = 3, additionally assume that (Hi) of Lemma \2.1\ 
holds. Replace v by UAr(-) = N'^^^v{ - N) and let (/i^"', . . . , /i^') be any min- 
imis er for the ground product state energy. Define (f)]^ — jj X^i^iC^i"')^- Then 
we have 

^lim^Xr=Xrr4 and - (C^;)'- 

where the convergence of 0^ is in the weak L^(U.'^)-sense and weakly for the 
probability measures (/^^^{x) dx towards the measure {(/>'- ^J^)'^ (x) dx . 

Note that, in d = 3, the interaction potential is rescaled in the same way in 
Theorems 12.31 and 12.41 However, the two relevant parameters depend on differ- 
ent properties of the potential (the scattering length, respectively the integral) 
and have different values, since a{v) < a{v) (see Section [^TT|) . In particular, 
for large enough, the ground state of xn is not a product state. This im- 
plies the strictness of the inequality for the two ground state energies, for v 
replaced by vn{-) = N^v{-N). The phenomenon that (unrestricted) ground 
states are linked with the scattering le ngth has been theoretically predicted 
for more genera l TV-body problems (see Fetter. A.L. and Walecka. J.D.I ( 197lL 



Indeed, Landau combined a diagrammatic method (a 



Born approximation of the scattering length) with Bogoliubov's approximations 
to almost reconstruct the scattering length from the L^-norm of w o | ■ | in the 
(non-dilute) ground state. However, the relation between the L^-norm and the 
product ground states was not rigorously known before. 

In d = 2, a more substantial difference between the large- A^ behaviours of 
Xn and Xjv' is apparent. Not only the asymptotic relation between the reach 
and the strength of the repulsion is different, but also the order of this rescaling 
in dependence on N . We can offer no intuitive explanation for this. 

Interestingly, in the hard-core case, Xjv' shows a rather different large-A^ 
behaviour, which we want to roughly indicate in a special case. Assume that W 
and V are purely hard-core potentials, for definiteness we take W = cxd115j(o)<= 
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and V = oo]l[o_a]- We replace v by vm{-) — v{ ■ /^n) for some ^at i (a pre-factor 
plays no role). Then Xn'' equal to times the minimum over the sum of 
the principal Dirichlet eigenvalues of —A in subsets of the unit ball having 
distance > a/3jv to each other, where the minimum is taken over the A'^ sets. It 
is clear that the volumes of these N sets should be of order , independently 
of the choice of ^n- Then their eigenvalues are at least of order TV^/'^. Hence, 
one arrives at the statement liminfjv->oo ^^^^'^Xjv' Xn^ tends to oo 

at least like iV^/'^. 

3 The probabilistic models 

Much thermodynamic information about the Boson system is contained in the 
traces of the Boltzmann factor e"''^" for (3 > 0, like the free energy, or the 
pressure. Since the 1960ies, interacting Brownian motions are generally used 
for probabilistic representations for these traces. The parameter which is 
interpreted as the inverse temperature of the system, is then the length of the 
time interval of the Brownian motions. 

However, the traces do not contain much information about the ground state. 
Since the pioneering work of Donsker and Varadhan in the early 1970ies it is 
basically known that the ground states are intimately linked with the Brownian 
occupation measures. This link is rigorously established via the theory of large 
deviations for diverging time, which corresponds to vanishing temperature. 
We introduce two different models of interacting Brownian motions. These 
models are given in terms of transformed measures for paths of length /3 in terms 
of certain Hamiltonians. Let a family of N independent Brownian motions, 
{Bl^^)t>o, . . . , (i?f"')f>o, in M"* with generator —A be given. The Hamiltonians 
of both models possess a trap part and a pair- interaction part. The trap part 
is for both models the same, namely 

Hn..i3 = Y1 / W{Bi^^)ds. (8) 

The Hamiltonian of our first model consists of two parts: the trap part given in 
([5]), and a pair-interaction part, 

l<i<j<N 

We look at the distribution of the N Brownian motions under the trans- 
formed path measure 

dPjv,^ = -J— exp{-HN.p-GN,i3)dP, where Z^./j = ¥.(exp{~HN^i3-GN,i3)) ■ 

Here E denotes the Brownian expectation for deterministic start at the origin 
and time horizon [0,/3]. We call Fn,p the canonical ensemble model, since it 
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is derived, via a Feynman-Kac formula, from the trace-class operator of the 
canonical ensemble, e~^^" . That is, the trace is given as 



Tr(e 



da;i 



N 



da; TV I 



'NJ3 



Here denotes the expectation with respect to a Brownian bridge starting 

in Xi and terminating in Xi after time 

However, a system of N Bosons is described by a trace of the projection to 
symmetric wave functions, i.e., wave functions that are invariant under permu- 
tations of the single particle indices. Hence the trace for a system of Bosons is 
given as 



Tr + (e 



) = - y 



dxi 



\ Ax^^^l U-"«-^- 

^=l ^ 



where Sjv is the group of permut ation of TV elein ents. These symmetrised sys- 
tems are the subject o f the review Adams! (l2007bl) in these proceedings. Recent 
result s can be found in Adains. S. and Dorlas. T.l ( 2007a ). Adams. S. and Konig. W 
(2003) andlAdamd (l2007d) . lAdamd (l2007dl) . 

The path measure Pat,/? is a model for N Brownian motions in a trap W with 
the presence of a repellent pair interaction. We can conceive the A^-tuple of the 
motions, Bt = (B, 



(1) 
t ) 



Bl ), as one Brownian motion in 



normalised occupation measure of the dA^-dimensional motion, 

/3 



1 r 

^p{dx) = - I 5B,(dx)ds, 
Jo 



and 



K 



I) dsAt. 



(9) 



E ^ / / Ki^, 

Note that the i-th Brownian motion interacts with the mean of the whole 
path of the j-th motion, taken over all times before (i. Hence, the interac- 
tion is not a particle interaction, but a path interaction. The interaction ([9|) 



Introduce the 



which is a random element of the set A^i(R''^) of probability measures on 
W^^ . It measures the time spent by the tuple of N Brownian motions in a given 
region. Note that there is only one time scale involved for all the motions, i.e., 
the Brownian particles interact with each other at common time units. We can 
write the Hamiltonians in terms of the occupation measure as 



Note that the energy functional (h^Tij^h) may be rewritten {h,Ti.]^h) = 
/jv(^) -I- (2U, /i) + (d, fi) for the probability measure fi{dx) — h'^{x) dx. 
Our second Brownian model is defined in terms of another Hamiltonian. We 
keep the trap Hamiltonian H^^fj as in ([5]), but the interaction Hamiltonian is 
now 
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is related to Polaron type m odels Donsker. M.D. and Varadhan. S.R.S. ( 1983f ). 
Bolthausen. E. et alj (|l993r ). where instead of several paths a single path is con- 
sidered. We consider the corresponding transformed path measure, 



dP5^,';3 = ^(«r exp(-Hjv,/3-ifjv,/3) dP, where Z^j^ = E{exp{-HNj3-KNj3)) . 

In Theorem 14.31 below it turns out that the large-/? behaviour of Z''^''^ is inti- 
mately related to the Hartree formula in fT]). Therefore, we call this model the 
Hartree model. At the end of this section we comment on its physical relevance. 
We introduce the normalised occupation measure of the i-th motion, 

1 

P Jo 

The tuple of the N occupation measures, (/i^^', . . . , /^Jj"'); plays a particular role 
in this model. We can write the Hamiltonians as 

Fw,/3=/3(2U,m|) and Kn,p = f3 ^ ,V^lf) ^ tif), 

l<i<j<N 

where we recall the operator V with kernel w o [ • |, and /i^ = ^jj^' (JD ■ • • ® /x^"' is 
the product measure. 

4 Large deviations results 

We present our main large deviations results for both the canonical ensemble 
and the Hartree model. In Section |4T] the zero temperature (i.e., large-/?) limit 
is considered, and in Section 14.21 the large-A^ limit, both at zero temperature 
and positive temperature. 

4.1 Vanishing Temperature 

It turns out that the large-/? behaviour of the canonical ensemble model is de- 
scribed by the ground state of the Hamilton operator Tijv via a large deviations 
principle. The rate function /jy appearing in Theorem 14.11 is the well-known 
Donsker- Varadhan rate function on M''^ defined by 

/.(.) = I ll^Vltli; if vlfei^^(M-) exists, ^^^^ 
[ oo otherwise. 

Simplifying, the large deviations principle says that, as /3 — > cxd, 

P(^^ «^) «e-^^"(^\ neMiiR''^). 

Theorem 4.1 (Canonical ensemble model at late times) Fix N G N. 
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ir7^1ogIE(exp(-ffAr,^ - Gn.p))= -Xn, 

where xn is the ground-state energy per particle of the N -particle operator 
Hn given in ([6]). 

(a) As (3 — > oo, the distribution of up on A4i(R'^^) under Pat,/? satisfies a 
principle of large deviation with speed (3 and rate function In given by 

/^(m) = /Ar(M) + (22J,/i) + -NxN for^i G MiiR''''). 

(Hi) The distribution of fip under P^^p converges weakly towards the measure 
h^,{x)'^ dx, where h^, is the unique minimiser in 



Remark 4.2 It is well-known I Ginibrt 1 197(\ ) that the bottom of the spectrum 
ofUM is related to the large-j3 behaviour of the trace o/e^'^^", more precisely, 

Xn ^ - lim — logTr (e"'''""). 

O 

Theorem 4.3 (Hartree model at late times) Assume that W and v are 
continuous in {W < oo} resp. in {v < oo}. Furthermore, assume in the soft- 
core case that there exists an e > and a decreasing function v: (0, e) — > M with 
V <v on (0, e), which satisfies Jg j.^^ G(0, y)v{\y\) dy < oo, where G denotes the 
Green's function of the free Brownian motion on . Fix iV G N. 

lim J^logE(exp(-i?jv,/3 - K^.p))^ -xT ■ 

(ii) As P oo, the distribution of the tuple {Hp\ ■ ■ ■ of Brownian oc- 

cupation measures on A^i(R'')^ under P^^'^ satisfies a large deviation 
principle with speed (3 and rate function 

N 
i=l 

with /ii, . . . ,fj,N G Aii{M.'^) where Ii is defined in (fTU|) . and jj.® = ® 
■ • • ® is the product measure. 

(Hi) The distribution of {fi'j^\ . . . , fi'j^') under P^^'^ is attracted by the set of 
minimisers for ground product state energy . 
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4.2 Large systems at Positive Temperature 

We now formulate our results on the behaviour of the Hartree model in the limit 
as N ^ oo, with f3 > fixed. As in the zero temperature case in Theorem l2.41 we 
replace v by V]\{{-) = N'^~^v{-N); we write -ft^j^^ for the Hamiltonian introduced 
in ®. 

First we introduce an important functional, which will play the role of a 
probabihstic energy functional. Define Jp: Mi{'R.'^) — + [0, oo] as the Legendre- 
Fenchel transform of the map Cb(R'^) 3 / ^ ilogE[c/o^(^=)d"] on the set 
Cb(K'') of continuous bounded functions on M"*, where {Bs)s>a is one of the 
above Brownian motions. That is, 

sup f(/i,/)-ilogE(e^o'-^(^=)<^^)}, fieMiiR"). 

/GCb(R<*) P ' 

Here A^i(R'') denotes the set of probability measures on . Note that J/? 
depends on the initial distribution of the Brownian motion. One can show that 
J/3 is not identical to +oo. Clearly, is a lower semi continuous and convex 
functional on A4i(K''), which we endow with the topology of weak convergence 
induced by test integrals against continuous bounded functions. However, 
is not a quadratic form coming from any linear operator. We wrote (/i, /) = 
/ifd /(s;) /i(da;) and use also the notation (/, g) = J^^ f{x)g{x) dx for integrable 
functions /, g. If ^ possesses a Lebesgue density c^^ for some L^-normalised 
(j) e L"^, then we also write Jpi^fj)^) instea d of Jaitx). It turns out th at J/3(^) — oo 
if /i fails to have a Lebesgue density, see lAdams. S. et al. ] (|2006bl) . 



In the language of the theory of large deviations, Jp is the rate function that 
governs a large deviations principle. The object that satisfies this principle is 
the mean of the N normalised occupation measures, 



1 ^ 



1=1 



Roughly speaking, this principle says that, as ^ oo, 

The principle follows from Cramer's theorem, together with the exponential 
tightness of the sequence (M7v,/3)A'eN- 

To apply this principle, we have to express our Hamiltonians H^^p and K^^p 
as functionals of /Iat g- For the first this is an easy task and can be done for any 
fixed A^: 

/ I ^ 

Now we rewrite the second Hamiltonian, which will need Brownian intersec- 
tion local times and an approximation for large A". Let us first introduce the 



r 1 

Nf3 / Wix)-Y,f^p\<i^) = Np{W,-p^^p). 
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intersection local times, see lGeman. D. et al. (1984). For the following, we have 
to restrict to the case d G {2, 3}. 

Fix 1 < i < j < N and consider the process B^'^—B^'\ the so-called confluent 
Brownian motion of _B''' and —B'^'K This two-parameter process possesses a 
local time process, i.e., there is a random process {L'-^''\x))^fzRd such that, for 
any bounded and measurable function / : R'' R, 

/ f{x)L^^-\x)dx f ds f dtf{Bf - B'f) 

JB<i P Jo Jo 

ji^'p'^\dx)4'^\dy)f{x-y). 



Hence, we may rewrite K^^''^ as follows: 

l<i<i<N 



It is known (jGeman. D. et all . 11984 Th. 1) that i^''' {x))^^Ts.d may be chosen 
continuously in the space variable. Furthermore, the random variable L^''^'(0) = 
liuix^o L'j^''\x) is equal to the normalised total intersection local time of the 
two motions B^'' and up to time /3. Formally, 



ds 



]1{B« e dx} 







dx 



/3 



dt 



t{Bi'^ e da:} 
dx 



(i) 

U 

dx 

A dx 



dx 



Using the continuity of Lj^' , we approximate 



N N 



l<i<j<N 



dx 



Hence, using Varadhan's lemma and ignoring the missing continuity of the 
map /i I— !■ Il^llii this heuristic explanation is finished by 

:(cxp { - Nf3[{W - f,-pN,p) ~ ^Mv) 



E 



dx 



}) 
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where 

Xr(/3)= , inf„ „ {Mc^') + {W,^')+'inaM\i}. (11) 

0GL2(Rd): ||0||2 = 1 I- J 

Here we substituted (/)^(a;)dx — /i(da:), we may restrict the infimum over 
probabihty measures to the set of their Lebesgue densities 0^ . 
Let us now give the precise formulation of our results. 

Theorem 4.4 (Many-particle limit for the Hartree model) Assume that 
d e {2,3} and let W and v satisfy Assumptions (W) and (v), respectively. In- 
troduce 

a{v) :— I v{\y\)Ay<oo. 

Fix f3 > 0. Then, as N ^ oo, the mean Jl^ p = jj 12iLi Z^/?' '^f normalised 
occupation measures satisfies a large deviation principle on Ali(R'*) under the 

-H -if'"' 

measure with density e with speed N (3 and rate function 

= [Jpi^^) + (W^'-^') + ll<^lll ^f4>^^^ exists, 

^ 1 00 otherwise. 

The level sets 6 Mi{W^): /^®'(Ai) < c}, c € M, are compact. 
Lemma 4.5 (Analysis of X^'(/3)) Fix /3 > and a > 0. 

(i) There exists a unique Lp' -normalised minimiser 0* g L/^iW^) n L^{W^) of 
the right hand side of (llip . 



(a) For any neighbourhood N C L'^{W^) n L'^{W^) of(j)^, 

inf {j^(0^) + (W^,02}+47ra||0|||}>xr(/3). 

Here 'neighbourhood' refers to any of the three following topologies: weakly 
in , weakly in L'^ , and weakly in the sense of probability measures, if (p 
is identified with the measure (f>{x)'^ dx. 

Corollary 4.6 (Free energy for positive temperature) Let the assumptions 
of the previous Theorem \4.4\ be satisfied. Then the specific free energy per par- 
ticle is 

lim logE(e"-^"''^"-^""3) =y<®' (/?). 

References 

Adams, S. (2007a). Interacting Brownian bridges and probabilistic interpreta- 
tion of Bose-Einstein condensation. Habilitation thesis, University Leipzig. 



17 



Adams, S. (2007b). Large deviations for empirical measures in cycles of integer 
partitions and their relation to boson symmetrisation. 

Adams, S. (2007c). Large deviations for empirical path measures in cycles of 
integer partitions, preprint arXiV:math.PR/0702052v2. 

Adams, S. (2007d). Large deviations for shape measures of integer partitions 
under non-uniform distribution, in preparation. 

Adams, S. and Bru, J.-B. (2004a). Critical Analysis of the Bogoliubov Theory 
of Superfluidity. Physica A 332, 60-78. 

Adams, S. and Bru, J.-B. (2004b). Exact solution of the AVZ-Hamiltonian in 
the grand-canonical ensemble. Ann. Henri Poincare 5, 405-434. 

Adams, S. and Bru, J.-B. (2004c). A New Microscopic Theory of Superfluidity 
at all Temperatures. Ann. Henri Poincare 5, 435 476. 

Adams, S., Bru, J.-B., and Konig, W. (2006a). Large deviations for trapped 
interacting Brownian particles and paths. Ann. Probab. 5^(4), 1340-1422. 

Adams, S., Bru, J.-B., and Konig, W. (2006b). Large systems of path-repellent 
Brownian motions in a trap at positive temperature. Electronic Journal of 

Probability 11. 

Adams, S., CoUevecchio, A., and Konig, W. (2007). A variational formula for 
the free energy of non-dilute many-particle systems, in preparation. 

Adams, S. and Dorlas, T. (2007a). Asymptotic Feynman-Kac formulae for large 
symmetrised systems of random walks. Annales de I'institut Henri Poincare 
(B) Probabilites et Statistiques . in press, preprint arXiV:math-ph/0610026. 

Adams, S. and Dorlas, T. (2007b). C*-algebraic approach to the Bose-Hubbard 
model, preprint math-ph arXiv: 0708.0518v. 

Adams, S., Georgii, H.O., and Kotecky, R. (2007). Gibbs measures for random 
permutations, in preparation. 

Adams, S. and Konig, W. (2007). Large deviations for many Brownian bridges 
with symmetrised initial- terminal condition. Prob. Theory Relat. Fields, ac- 
cepted for publication, preprint arXiV:math. PR/0603702. 

Anderson, M.H., Ensher, J.R., Matthews, M.R., Wicman, C.E., and Cornell, 
E.A. (1995). Observation of Bose-Einstein condensation in a dilute atomic 
vapor. Science 269, 198-201. 

Bolthausen, E., Dcuschcl, J.D., and Schmock, U. (1993). Convergence of path 
measures arising from a mean field or polaron type interaction. Probab. Theory 
Related Fields 95, 283-310. 

Bose, S. (1924). Zeitschrift fiir Physik 26, 178. 



18 



Bradley, C.C., Sackctt, C.A., Toilet, J.J., and Hulct, R.G. (1995). Evidence 
of Bose-Einstein condensation in an atomic gas with attractive interactions. 
Phys. Rev. Lett. 75, 1687-1690. 

Dalfovo, F., Giorgini, S., Pitaevskii, L.P., and Stringari, S. (1999). Theory of 
Bose-Einstein condensation in trapped gases. Rev. Mod. Phys. 71:3, 463-512. 

Davis, K.B., Mewes, M.O., Andrews, M.R., van Drutcn, N.J., Durfee, D.S., 
Kurn, D.M., and Kcttcrlc, W. (1995). Bosc-Einstcin condensation in a gas of 
sodium atoms. Phys. Rev. Lett. 75(22), 3969-3973. 

Dickhoff, W.H. and Van Neck, D. (2005). Many-Body Theory Exposed. World 
Scientific, Singapore. 

Donsker, M.D. and Varadhan, S.R.S. (1983). Asymptotics for the polaron. 
Comm. Pure Appl. Math. 36:4, 505-528. 

Einstein, A. (1925). Sitzher. Kgl. Preuss. Akad. Wiss. 3. 

Fetter, A.L. and Walecka, J.D. (1971). Quantum Theory of Many Particle 
Systems. McGraw-Hill, New York. 

Feynman, R. (1953). Atomic theory of the A transition in Helium. Phys. Rev. 91, 
1291-1301. 

Feynman, R. (1955). Progress in Low Temperature Physics, Vol. L, pp. 17. 
Niorth Holland Amsterdam. 

Geman, D., Horowitz, J., and Rosen, J. (1984). A local time analysis of inter- 
section of Brownian paths in the plane. Ann. Probab. 12:1, 86-1070. 

Ginibre, J. (1970). Some Applications of Functional Integration in Statistical 
Mechanics and Field Theory. Gordon and Breach Science Publ., New York. 

Gross, E. (1961). Structure of a quantized vortex in boson systems. Nuovo 
Cimento 20, 454-477. 

Gross, E. (1963). Hydrodynamics of a superfluid condensate. J. Math. Phys. 4:2, 
195-207. 

Landau, L. (1941). The theory of superfluidity of Helium II. J. Phys. USSR 5, 
71. 

Landau, L.D. and Lifshitz, E.M. (1951). Statisticheskai Fizika. Fizatgiz Moscow. 

Lieb, E.H. and Seiringer, R. (2002). Proof of Bose-Einstein condensation for 
dilute trapped gases. Phys. Rev. Lett. 88:17, 170409-1-4. 

Lieb, E.H., Seiringer, R., and Yngvason, J. (2000a). Bosons in a trap: A rigorous 
derivation of the Gross-Pitaevskii energy functional. Phys. Rev. A 61, 043602- 
1-13. 



19 



Lieb, E.H., Sciringcr, R., and Yngvason, J. (2000b). The ground state energy 
and density of interacting Bosons in a trap, in: Quantum Theory and Symme- 
tries, Goslar 1999; H.D. Doebner, V.K. Dobrev, J.D. Hennig and W. Luecke, 
eds., World Scientific, 101-110. 

Lieb, E.H., Seiringer, R., and Yngvason, J. (2001). A rigorous derivation 
of the Gross-Pitaevskii energy functional for a two-dimensional Bose gas. 
Comm. Math. Phys. 221 17-31. 

Lieb, E.H. and Yngvason, J. (2001). The ground state energy of a dilute two- 
dimensional Bose gas. J. Statist. Phys. 103, 509. 

London, F. (1938). The A-Phenomenon of Liquid Helium and the Bose-Einstein 
Degeneracy. Nature 14 1, 643. 

Onsagcr, L. (1949). Statistical hydrodynamics. Nuovo Cimento 6(Suppl. 2), 
279-287. 

Onsager, L. and Penrose, O. (1956). Bose-Einstein condensation and liquid 
Helium. Phys. Rev. 101 576-584. 

Penrose, O. (1951). On the quantum mechanics of Helium II. Phil. Mag. ^2, 
1373-1377. 

Pitaevskii, L. (1961). Vortex lines in an imperfect Bose gas. Sov. Phys. JETP 13, 
451-454. 

Pitaevskii, L. and Stringari, S. (2003). Bose-Einstein Condensation. Clarendon 
Press, Oxford. 

Popov, V. (1983). Functional Integrals in Quantum Field Theory and Statistical 
Physics. Ricdcl, Dordrecht. 

Ruelle, D. (1969). Statistical Mechanics: Rigorous Results. Addison- Wesley. 



20 



